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I. INTRODUCTION 

The discovery of neutrino oscillations triggered a nat- 
ural interest in the problem of masses of neutrinos. Un- 
fortunately, in the oscillation experiments only the dif- 
ferences of the squares of masses can be determined. 
The importance of knowing the absolute scale of neu- 
trino masses is obvious. This knowledge will allow to 
set direction in which the changes of the standard model 
of particles and interactions should go; what is more, 
the problems of dark matter and dark energy, lepto- and 
baryogenesis, evolution of the Universe and many other 
could be addressed. Even if the values of neutrino masses 
will be measured in experiments, the question about the 
mechanism of obtaining those masses remains open, just 
in the same way as the widely approved Higgs mechanism 
is still not experimentally confirmed. 

There are many proposals of generating the neutrino 
mass matrix. Starting from ad-hoc ansatze, through the 
most-widely approved see-saw model, through extra di- 
mensions, through a result of supersymmetry breaking. 
Among these one also finds the loop mechanism for Ma- 
jorana neutrinos. The effective vertex of the form Diy 
is there expanded to contain a squark-quark or slepton- 
lepton loop. This setting introduces the R-parity viola- 
tion, therefore it needs to be described within a super- 
symmetric model with explicitly or spontaneously broken 
R-parity. Such models provide an elegant way of not only 
resolving the naturally small neutrino mass problem, but 
also introduce supersymmetry, needed by the string the- 
ory, solve the hierarchy problem, provide much better 
description of the anomalous magnetic moment of the 
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muon and much more. 

The problem of R-parity violation in supersymmet- 
ric models received a great deal of attention during last 
few years The studies of this topic were connected 
with, among others, leptonic decays gravitino decays 
y| and the problem of neutrino masses and oscillations 
in IS' In general, supersymmetric models with R- 
parity violation (RpV) fall into one of the three cate- 
gories. First, we have the spontaneous RpV. In this case 
the R-parity is violated by a non-zero vacuum expecta- 
tion value of some scalar field 0. Another possibility 
is the explicit breaking of R-parity by introduction of 
bilinear and/or trilinear terms. The bilinear RpV mod- 
els 8] are characterized by good predictivity due to a 
small number of parameters. The third category is the 
explicit RpV by trilinear terms present in the superpo- 
tential j^, 0| . In this case one allows for the presence of 
bilinear terms, since these would anyway show up dur- 
ing the RGE evolution of trilinear coupling constants, 
assuming at the same time that they do not affect the 
phenomenology of the trilinear terms. This is motivated 
by the fact that there is no fine tuning among differ- 
ent contributions, which can therefore be analyzed sepa- 
rately. The trilinear scenarios are the most studied due to 
the reachest phenomenology and possibility of obtaining 
most interesting limits on non-standard physics parame- 
ters. 

Our work follows the line of research concerning the 
generation of neutrino masses in SUSY without R-parity, 
that has been developed in the last few years. The aim of 
our paper is to get new individual limits on the R-parity 
breaking parameters by taking the advantage of the re- 
cent data on neutrino oscillations and neutrinoless double 
beta decay. In addition, the previous studies presented 
in |y, lid III ] are improved by a more accurate treatment 
of the neutrino mass contributions, in particular by re- 
ducing the dependence on the SUSY parameter space. 
The considered model is the minimal supersymmetric 
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standard model with supersymmetry breaking transmit- 
ted by (super)gravity interactions (SUGRA MSSM) 
with the squark and slepton mixing phenomena properly 
included. The RGE evolution within GUT constrained 
SUGRA MSSM is introduced to obtain the low-energy 
particle spectrum. The GUT constraints involve unifying 
masses and coupling constants to some common values 
at the GUT scale. These are universal scalar mass mo, 
gaugino mass rai/2, trilinear scalar coupling Aq, the ra- 
tio of the Higgs vacuum expectation values tan (3 and the 
sign of the bilinear Higgs mixing parameter /i. 

The paper is organized as follows. In the next sec- 
tion we present the model and describe our procedure 
of finding low energy spectrum of SUSY particles using 
the GUT constraints. We present also improved versions 
of contributions of different loops to the neutrino mass 
matrix. In Sec. Ill we present upper limits on various 
combinations of RpV coupling constants, both the trilin- 
ear A's and dimcnsionful bilinears A's. Discussion and 
conclusions follow at the end. 



II. THE MODEL 

The loop mass mechanism may be described in the 
framework of R-parity violating MSSM (RpVMSSM) 
with trilinear and bilinear soft breaking terms. The 
MSSM (see e.g. ^^'^ its many variations are well 

known in the literature. Here, we closely follow Ref. 
regarding the conventions and construction of the mass 
matrices. In short, RpVMSSM is characterized by the 
superpotential which consists of the R-parity conserving 
part 



W 



MSSM 



(1) 



and the R-parity violating part 



W^- = eab 
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-XijkL^LjEk + \ijf^L°lQj Dkx 



+ ^exyAl.kUrDpl + eabK^LfHl (2) 

The Y's are 3x3 Yukawa matrices. L and Q are the 
SU{2) left-handed doublets while E, U and D denote the 
right-handed lepton, up-quark and down-quark SU (2) 
singlets, respectively. Hi and H2 mean two Higgs dou- 
blets. We have introduced color indices x,y,z — 1,2,3, 
generation indices i, j,k — 1,2,3 and the SU(2) spinor in- 
dices a,b,c = 1,2. In order to get rid of too rapid proton 
decay and to describe lepton number violating processes, 
like the neutrinoless double beta decay, it is customary 
to set A" = 0. 

We supply the model with scalar mass term 

j^mass ^ mjj^h\hi + mjj^hlh2+q^mlq + l^mll 

+ wm^ut _^ dmjyd^ + emle^, (3) 



soft gauginos mass term 

^gaug. ^ 1 (^MiB^ B + N^W^ + Msgjg" + h.c. 

'(4) 

as well as the supergravity mechanism of supersymmetry 
breaking, by introducing the Lagrangian 
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{Au),,qrhlu,x + Btihihl + B2ed'ih% (5) 



where lowercase letters stand for scalar components of 
respective chiral superfields, and 3x3 matrices A as well 
as B^ an B2 are the soft breaking coupling constants. 

All the running parameters are obtained by using the 
renormalization group equations (RGE) jlj, [T^ . At the 
beginning, one evolves all gauge and Yukawa couplings 
for three generations up to the GUT scale Mgut ~ 
10^^ GeV. We use the one-loop standard model RGE [l^ 
below the mass threshold, where SUSY particles start to 
contribute, and the MSSM RGE jl3| above that scale. 
The contribution of 2-loop diagrams as well as those com- 
ing from the R-parity violating couplings has been proven 
to be irrelevant in discussions such as ours The 
SUSY scale is initially set to 1 TeV for all particles and 
is dynamically modified together with evolution of their 
masses. At the GUT scale the masses of all scalars and 
fermions are set to a common value mo = mi/2 = tu- We 
have considered a "small", chosen to be 150 GeV, and a 
"large" (1000 GeV) value of m in our analysis. We unify 
also the soft trilinear couplings according to — AoYi, 
with Aq ~ 500 GcV. We postpone the discussion of the 
influence of ^0 on the results to a forthcoming paper. In 
the next step we construct all the relevant mass matrices 
(squark, slepton, chargino and neutralino) and perform 
RGE evolution of all the quantities back to Mz scale, 
taking care of the minimization of the tree-level Higgs 
potential (important for EWSB breaking) and radiative 
corrections. After iterating this procedure and obtaining 
stable values of the parameters, we confront the obtained 
values with restrictions coming from the present theoret- 
ical assumptions and phenomenological data. Those con- 
straints involve (1) finite values of Yukawa couplings at 
the GUT scale; (2) proper treatment of electroweak sym- 
metry breaking; (3) requirement of physically acceptable 
mass eigenvalues at low energies; (4) FCNC phenomenol- 
ogy. 

The first problem is related to the values of tan f3 and 
is checked during the RGE procedure. For very small 
tan/3 (< 1.8) the top Yukawa coupling may "explode" 
before reaching the GUT scale. It follows from the fact 
that Ytop{Mz) ^ 1/sin/?. Similarly, other couplings Yi, 
and Yr "blow up" before the GUT scale for tan^ > 50 
because they are proportional to 1/ cos f3 at electroweak 
scale. In our analysis we have kept tan/3 « 20 leaving 
the detailed discussion to a forthcoming paper. 

Another theoretical constraint is imposed by the 
EWSB mechanism. In order to obtain a stable min- 
imum of the scalar potential, the following conditions 
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must hold: 

2Bfi < 2\^if +m%^^+ni%^. (6) 

They are always checked in our procedure during RGE 
running, and points which do not fulfill these conditions 
are rejected. Next restriction comes from the require- 
ment of positive eigenvalues of mass matrices squared at 
the electroweak scale. The last requirement (see, e.g., 
|l5j for details) comes from the strongly experimentally 
suppressed FCNC processes and provides the most severe 
constraints. 

The so-obtained low energy spectrum is then used in 
further calculations. Although, as will be seen, only the 
squarks and sleptons masses enter the formulas, they de- 
pend in a complicated way on all other masses and cou- 
pling constants through RGE equations ^V?\. Therefore 
a complete and careful treatment is necessary. 

The neutrino mass matrix consists in our approach of 
three main parts: 

M^ = M''''"' + M' +M'^, (7) 

which are the tree level value and the contributions com- 
ing from lepton-slepton and quark-squark loops, respec- 
tively. We note that there are other terms that may be 
included in Eq. Q, in particular loops contributions 
with bilinear insertions For the sake of simplicity, 
however, we do not consider them by following the popu- 
lar approach to get individual limits on R-parity breaking 
parameters. We adopt the conventional hypothesis that 
different contributions do not significantly compensate 
each other and for this reason it is possible to extract 
limits on individual contributions without knowing the 
others. 

Let us first recall the well known results. In the lowest 
order, the contribution to the mass matrix reads |10j 

Ml + M2 tan^ 0w 
^ AiliM^iMi + M2 tan2 Ow) sin 2/3 - MiAh^i^) ' 

(8) 

where = — {Hi)Ki, and {i>i) are the vacuum 

expectation values of the sneutrino fields. 

Beyond the tree-level, the Majorana neutrino mass ma- 
trix may also be generated by considering one loop self- 
energy diagrams. The particles that propagate inside the 
loops are either quark and squark or lepton and slepton. 
Let us start with the squark-quark loop. The relevant 
Feynman diagrams are shown on Fig. ^ It is important 
to note, that one may (and should, if one wants to be 
very accurate) consider not only the trilinear couplings, 
but also the mass insertions described by the bilinear 
terms in the superpotential and Lagrangian (see e.g. 0|). 
Here, however, we take the conventional approach and 
assume that the phenomenology of the trilinear contri- 
bution (which is more interesting from the point of view 
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FIG. 1: Feynman diagrams representing the squark-quark 
loop contribution to the Majorana neutrino mass. 
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FIG. 2: Feynman diagrams representing the slepton-Iepton 
loop contribution to the Majorana neutrino mass. 

of obtaining constraints on the RpV coupling constants) 
remains unaffected by the contribution coming from bi- 
linear terms. Contrary to earlier approaches 
we take into account the down-squark mixing exactly 

dh — di cos -t- 1^2 sin 6*, 

du, = —di sin 6* -I- ^2 cos 6*. (9) 

Here L and R label the left- and right-handed squark 
states in the weak basis, while the 1 and 2 subscripts 
denote the two mass eigenstates. The mixing angle is 
defined by 

sin(26''') = 2mqk{Ak + fJ,ta.n(3) 

X [(mL - mi -0.34Af|cos(2^))2 
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4to + ^tan/3)] 



(10) the three angles is 



with Ak = {Au)kk and tan/3 being the ratio of Higgs 
vacuum expectation values. The squark mass eigenstates 
take the forms 



- (m% + m% ) + m„j (m„j 



Aj + /X tan P 
sin(26'J) 



MzC0s(2/3), 



m% = - (m% + m^j ) + m„j f TO„j 



Aj + II tan /3 
sin(26'J) 



-M2C0s(2/3). 



(11) 



By introducing two dimensionless quantities, xj*" — 
m^, /m?fc and xi'^ = m}Jrn}-.k-, one arrives at the fol- 

' ^ q^ ' q^^ 

lowing form of the neutrino mass matrix 



——A' A' 



sin{2d'')m 



log(xf) , (xf + l)log(:rf) 



Jk 



I jk 



l)(xf -1), 



(12) 



We may repeat the same calculation for the slepton- 
lepton loop (see Fig. ISJ, just replacing in all definitions 
the squark masses and mixing by analogous quantities 
for sleptons, as well as the quark masses m^j with lepton 
masses m^j . The only difference will be the lack of the 
factor 3, which came in the previous case from summing 
over the three colors of quarks and, of course, different 
coupling constants. We end up with 



M'' , = 



1 



167r2 



^ijk^'kj sin(20'')mej 



]k 
V2 



i)fe 



k) 



(13) 



where now is the slepton mixing angle, j/^'^ = rn^j/m-f, 
and yl'' = mlJm?. 

'2 

Let us now explain the procedure for finding con- 
straints on the various products of coupling constants 
A, A', and A. The right hand sides of Eqs. (C2l, 
and (^21 can be calculated from the MSSM RGB run- 
nigs, during which the low energy particle/sparticle spec- 
trum is generated. We use random scatter to find sets 
of physically relevant values of the various parameters. 
In the next step the theoretical neutrino mass matrix 
is compared with the phenomenological three neutrino 
mass matrix in the flavor space, which is connected to 
the physical neutrino masses rrii by the mixing matrix U 
through the relation AiP'^ = U ■ diag (mi, 1112,^23) ■ . 
The standard parameterization of the Pontecorvo-Maki- 
Nakagawa-Sakata (PMNS) unitary matrix U in terms of 



C12C13 S12C13 siae 

-S12C23 - C12S23S13 C12C23 - S12S23S13 S23C13 
S12S23 - C12C23S13 -C12S23 - S12C23S13 C23C13 


















g»a2l 







^0 








(14) 



where Sy = sin 9ij , Cy = cos 9ij , and 9ij is the mixing an- 
gle between the flavor eigenstates labeled by indices i and 
j . The recent global analysis of neutrino oscillations 
yields the best fit values: sin^6'i2 = 0.3, sin^6'23 = 0.5 
and sin^6'i3 = 0.002. Note that for Majorana particles 
there appear three CP violating phases, one Dirac phase 
6 and two Majorana phases (q!2i and 031), which remain 
undetermined. Assuming the CP phases to be negligible 
one gets 



0.83 0.55 0.07 
U= I -0.42 0.55 0.72 
0.35 -0.63 0.69 



(15) 



The absolute scale of neutrino masses is not deter- 
mined by the neutrino oscillations, which depend only 
on differences of masses squared. From the global 
analysis [l9( of neutrino oscillations the best fit values 
Amii 6.9 10-^ eV^ and Am|i = 2.3 lO'^ eV^ are 
known. The three possible neutrino mass patterns are 
frequently considered |20|: 

i) The normal hierarchy (NH) of neutrino masses, which 
correspond to the case mi <^ m2 ^ TO3. Then we have 
nil ^ i/Amfj^, 77X2 — \/ Atoji and 7713 ~ ^yAm^. 

ii) Inverted hierarchy (IH) of neutrino masses. It is given 
by the condition <C mi < TO 2 . In the case for n eutrino 
masses we have 7713 <^ y^Am^ and mi ~ m2 ~ y'Arn^. 

iii) Almost degenerate neutrino mass spectrum: mi ~ 
m2 — ms. This case does not exclude the possibility 
that the lightest neutrino is much larger than -y/Aml^. 

The absolute scale of neutrino masses can be deter- 
mined by the observation of the end-point part of the 
electron spectrum of Tritium /3-decay, the observation of 
large-scale structures in the early universe and the detec- 
tion of the neutrinoless double beta decay (Oz^/3/3-decay) , 
if neutrinos are Majorana particles. The amplitude of 
the Oz//3/9-decay is proportional to the effective Majorana 
neutrino mass mpp — U^i mi + U^2 ™2 + U^3 "m-a- This 
process has not been seen experimentally until now and 
the best results have been achieved in the Heidelberg- 
Moscow (H-M) experiment (T^J^ > 1.9 x lO^^ y) 

(Recently, some authors of the H-M collaboration have 
claimed the experimental observation of the Oi^/3/?-decay 
of ^^Ge |23|. But the Moscow participants of the H- 
M collaboration, performing a separate analysis of the 
data, found no indication in favor of the evidence of the 
0j^/3/?-decay 0| . The disproof or the confirmation of the 
claim will come from future experiments.) By assuming 
the nuclear matrix element of Ref. |24| we end up with 
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|to/3/3| < 0.55 eV. With this additional input limit we 
can find the maximal allowed values for the matrix el- 
ements AA^^ of the neutrino mass matrix, which are as 
follows: 

/ 0.55 0.71 0.70 \ 

\ 0.70 0.70 0.76 J 

The elements of this matrix were obtained by assuming 
the whole allowed mass parameter space of neutrinos and 
all possible CP-phases of the neutrino mass eigenstates 
[Toj |. In the calculation we used the best-fit values of 
neutrino oscillation parameters given in Ref. |19| |. The 
elements of matrix (|16|l can be used to test various the- 
oretical approaches and allows one to extract limits on 
certain fundamental parameters. Of course, one can not 
expect that by diagonalizing of this matrix a relevant in- 
formation on the masses of neutrinos is obtained as each 
element of this matrix is a result of analysis of all pos- 
sible mixing of three neutrinos allowed by the neutrino 
oscillations and the 0i//3/3-deca.y data. 

Instead of taking into account the current limit on m^^ 
identified with the element A^fg we consider also other 
scenarios by assuming that the normal or inverted hier- 
archy of neutrino masses is realized in the nature. Then 
we get 



^j^ph~NH 


1 = 10-4 


eV 








/ (22.4 


- 27.2) 


(0.64- 


49.4) 


(5.16 


- 52.0) 


X (0.64 


- 49.4) 


(223 - 


273) 


(210 


- 267) 


\ (5.16 


- 52.0) 


(210- 


267) 


(196 


- 262) 




= 10-2 


eV 








( (1.86 


-4.72) 


(0.22- 


3.11) 


(0.26 


- 3.05) 


X (0.22 


-3.11) 


(0.62 - 


2.30) 


(0.96 


- 2.37) 


\ (0.26 


-3.05) 


(0.96 - 


2.37) 


(1.31 


- 2.50) 



These neutrino mass matrices were calculated by the as- 
sumption that the mass of the lightest neutrino is negli- 
gible (see the above definitions of the NH and the IH of 
neutrino masses). 

III. RESULTS AND CONCLUSIONS 

By confronting the phenomenological neutrino mass 
matrix M.p^^^^^ ^ derived from the analysis of the neu- 
trino data, with the theoretical mass matrix calculated 
within the R-parity breaking MSSM, it is possible to find 
constraints on various combinations of the lepton num- 
ber violating A, A' and A couplings, which enter Eqs. lO, 
(112), and linil- If M^''-^" and M^^^^" neutrino mass 
matrices are confronted with the theory, one ends up with 
predictions for the R-parity violation couplings. By con- 
sidering the maximal values of these matrices the largest 
possible values of R-parity breaking parameters are ob- 
tained. We note that the predictions for R-parity break- 
ing mechanisms associated with normal or inverted mass 



hierarchy are deduced by the assumption that one given 
mechanism dominates at a time. However, this scenario 
might be excluded by other phenomenology. We have 
used in our analysis the following quark masses: m„ = 5 
MeV, md = 9 MeV, = 175 MeV, = 1.5 GeV, 
rub = 5 GeV, rut = 174 GeV. 

Table ^ shows improved upper bounds on various 
combinations of coupling constants of the A, A' and A 
types, to be compared with the limits presented in Refs. 
[1, . The results were obtained for Aq = 500 GeV, 
mo = mi/2 = 150 GeV and 1000 GeV and for 
positive /i. We have found a weak dependence of the 
quantities under discussion on Aq SUSY parameter. Be- 
sides, we have kept tan/3 large (tan/3 « 20), leaving the 
discussion of the impact of this parameter on the results 
to the forthcoming paper. In general, the new bounds 
related to lower limit on the T^J^iJ^Ge) and neutrino os- 
cillation data are at least one order of magnitude stronger 
than those previously given It is mostly due to the 
assumption of the gravity-mediated (SUGRA) supersym- 
metry breaking and partially due to an improved treat- 
ment of the squark and slepton mixing. As expected the 
values of R-parity violating coupling related to normal 
hierarchy of neutrino masses are significantly suppressed 
in comparison with those related to the current lower 
limit on the Oz^/3/3-decay half- life . 

The new bounds are surprisingly close to those pub- 
lished in Ref. @, although the method used by the au- 
thors of these papers relayed on many simplifying as- 
sumptions. In particular, it involves setting some of the 
couplings to zero and assuming all other to be of the 
same order of magnitude. Also the whole mechanism of 
RGE running as well as GUT constraints were not used. 
In general the constraints in Q were Aj;33,A^33 < 10-® 
which is fully consistent with our results. The bounds on 
products of individual coupling constants in Tab. ^ are 
either of the same order of magnitude or 1-3 orders of 
magnitude stronger. 

A more optimistic scenario appears for the case of in- 
verted hierarchy. In general the corresponding values of 
product of A and A' coupling are about by factor four 
less stringent as those associated with the most strin- 
gent 0^/3/3-decay limit on the half-life. We stress again 
that these values of the R-parity breaking parameters 
were determined by the condition that a particular R- 
parity breaking mechanism dominates at a time. In some 
cases this might be excluded by the phenomenology of 
other processes. For example, from the R-parity break- 
ing SUSY mechanism of the 0:^/3/?-decay one gets the 
upper limit on the parameter X'm of the order of 10-^ 
|ll[|25l |. what is significantly less than the value presented 
in Table n 

In summary, we have used the GUT constrained R- 
parity violating minimal supersymmetric standard model 
to describe massive neutrinos. The three family neutrino 
mass matrix was calculated within framework including 
the tree-level neutrino-neutralino mixing and the one- 
loop radiative corrections. Then, the theoretical mass 
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TABLE I: Constraints on A, A' and A from their contribution to neutrino masses, using recent global analysis of the neutrino 
oscillation data [[l^, the currently best experimental limit on the Ou(3/3-decsLy half- life 21], and matrix element of Ref. |2^ . 
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matrix was compared with the phenomenological one, 
obtained by using the most recent global analysis of neu- 
trino oscillations data and the lower limit on the half-life 
of neutrinoless double beta decay of ''^Ge. This proce- 
dure allowed to improve the upper limits on certain prod- 



ucts of R-parity violating couplings, which are up to one 
order of rnaenitude more stringent as those previously 
published [Tl|. Further on, we assumed the normal and 
inverted hierarchy of neutrino masses and calculated the 
corresponding values of R-parity violating parameters of 
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the SUSY model under consideration. These values can 
be used in determining the perspectives of finding signal 
of R-parity violation in different experiments, in particu- 
lar at colliders. This issue is, however, beyond the scope 
of this paper. 
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